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A total Routhian surface (TRS) calculations were performed for the 
^*^Ce nucleus. The equilibrium deformation parameters and total angular 
momentum values have been calculated as a function of rotational fre- 
quency for both signatures of the lowest positive- and negative-parity ro- 
tational bands. Theoretical predictions are compared with the results of 
experimental studies of this nucleus. 

PACS numbers: 21.60.Ev, 27.60.-Hj, 21.10.Re 



The neutron-rich Ce nucleus vi^ith N=91 neutrons lies in a transitional 
region between nearly spherical lighter cerium isotopes and deformed ^^'^Ce 
with R=E4_|_/E2+ ~ 3.15 ratio[l] close to the rotational model value of 3.33. 
In N=91 nuclei Nilsson levels related to 217/2, II19/2 ^-ncl lii3/2 orbitals are 
located close to each other and to the Fermi surface at quadrupole defor- 
mation parameter of about /?2 ~ 0.2. Consequently, the rotational bands 
based on the above orbitals are expected to occur at low excitation ener- 
gies. Nilsson model calculations indicate, that the 2f7/2 and lhg/2 orbitals 
are strongly mixed in this deformation region. As a result of the mixing, 
the calculated decoupling parameters for 1/2~[541] state, a54i=0.65, and for 
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l/2~[530] state, a53o=0.15, are very different from values of about 4.6 and 
—3.4, respectively, expected for the pure states. With such low values of 
decoupling parameters the Coriolis mixing calculations performed with the 
use of the CORIOLIS code [2] show that all negative-parity rotational bands 
close to the Fermi level are normal strongly coupled bands with AI=1 se- 
quence and exhibit only small staggering of rotational levels. On the other 
side, for the positive-parity rotational bands close to the Fermi levels for 
N=91, the Coriolis mixing calculations predict existence of the decoupled 
rotational bands with AI=2 sequence of spins. 

To check these simple model predictions, an extensive total Routhian 
surface (TRS) calculations were performed for the ^^^Ce nucleus with the 
aim to compare theoretical predictions with the results of experimental stud- 
ies of the high spin states in this nucleus [3, 4] and the recent studies of the 
^^^La /3-decay carried out at the R2-0 reactor in Studsvik using the OSIRIS 
on-line fission product mass separator [5]. In the latter study tentative as- 
signment of spin and parity for several low-spin levels has been proposed 
and some ambiguities in the interpretation of the level scheme of ^^^Ce were 
resolved allowing for a comparison with theoretical predictions. 

A recently developed pairing-deformation self-consistent total Routhian 
surface (SC_TRS) model[6, 7] has been used. Within this model, the total 
Routhian of the nucleus (energy in the rotating frame) is calculated as a 
sum of macroscopic part and a microscopic correction accounting for shell 
effects and pair correlation. The liquid drop model of Ref . [8] and the shell 
correction method of Strutinsky[9] were employed. The total Routhian can 
be written as 

E'^(Z,N,/3) = E'^=0(Z,N,^) + [n^'iZ, N, P) - H''=^{Z, N, ^)} (1) 

where E'^^'^(Z,N,/3) represents the liquid drop energy, the single particle 
shell corrections, and pairing energy at frequency zero. The term in the 
brackets of Eq. (1) describes the change in energy due to rotation and $ 
denotes deformation parameters P2, Pa and 7. The total Routhian is cal- 
culated on a grid in deformation space, including quadrupole deformation 
(/32), triaxiality (7) and hexadecapole {f3i) degrees of freedom, and then 
minimized with respect to the shape parameters to obtain equilibrium pa- 
rameters. The deformed Woods-Saxon potential with an improved universal 
parameters set was used as a single-particle potential. In addition to the 
standard self-consistency requirement for the shape degrees of freedom, in 
this model the pairing correlations are also treated self-consistently. For 
more detailed description of the used TRS model see Ref. [10]. 

The calculated equilibrium deformation parameters of the positive-parity 
rotational band at fiuj ^ 0.075 MeV lie around the values of P2 = 0.228, 



3 



30 - 



HKi 20 - 



10 



■ 


exp 


(7C 




+ 


a = 


+ 1/2} 


• 


exp 






+ 


a ^ 


-1/2) 


▲ 


exp 


(JU 






a = 


-1/2} 




71 ^ 


+ , 


a 




+ 1/2 






K ^ 


+ , 


a 




-1/2 






n ^ 




a 




+ 1/2 






Jt = 




a 




-1/2 





0,0 



0,1 



0,2 0,3 

ho) [MeV] 



1 <- 



149 



Ce 



0,4 



0,5 



Fig. 1. The calculated and experimental values of the total angular momentum as 
a function of rotational frequency for both signatures of the lowest positive- and 
negative-parity rotational bands in ^"'^Ce. 



7 = -0.2° and = 0.085 while for the negative-parity rotational band 
these parameters are close to the values of (32 = 0.224, 7 = 0.1° and P4 = 
0.080. Both total Routhian surfaces minima are well defined and show no 
indication of 7-softness. They remain pretty stable with increasing rota- 
tional frequency up to Hlo a; 0.5 MeV. The main components of the wave 
functions of the lowest positive and negative parity bands are, respectively, 
the 3/2"'" [651] and 3/2~[521] configurations, in agreement with the Coriolis 
mixing calculations [5]. Relatively large calculated values of the hexade- 
capole deformation parameter, especially for the positive-parity band, and, 
simultaneously, low values of the 7 non-axial deformation are worth to note. 

In order to check the possibility that the obtained large values of the 
hexadecapole deformation parameter are generated due to neglecting higher 
order multipole expansion terms, e.g. (Sq, in our TRS model calculations, 
an additional calculations were performed minimizing potential energy of 
the lowest positive and negative states as a function of P2, P3, Pa, Pb and Pe 
deformation parameters at ^o; = MeV. The resulting P2 and P4 parameters 
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did not change significantly in the presence of liighcr multipole terms. Other 
calculated deformation parameters were very close or equal to zero, showing 
in particular that at low values of rotational frequency the octupolc degree 
of freedom plays negligible role for the lowest states in the ^"^^Ce nucleus. 
Only the 5/2"*" [642] configuration, which becomes dominant in the positive 
parity rotational band wave function at higher rotational frequencies {hoj > 
0.3 MeV), shows significant octupole deformation of P3 ^ 0.10 with 
0.03. 

The calculated values of the angular momentum <I-i;>as a function of 
rotational frequency Hlu are shown in Fig.l for both positive- and negative- 
parity rotational bands. Similarly to simple Coriolis mixing calculations a 
sizeable signature splitting is predicted only for the positive parity band. 
The backbending for the negative parity band is predicted at much lower 
frequency of about 0.20 MeV than for the positive parity band (about 0.33 
MeV). 

The experimental data (full circles and squares) for the rotational band 
built on the 3/2"*", 133.5 keV level with a new spin and parity assignments 
proposed in Ref . [5] match very well theoretical predictions only if the high 
spin structure represents the positive parity band. Consequently, it ex- 
cludes the possibility that the rotational band observed in experiment has 
a negative parity. Also the backbending frequency observed in the experi- 
ment (0.30 MeV), agrees only with prediction for the positive-parity band, 
although is about 0.03 MeV lower than the calculated value. Moreover, 
the only experimental point (full triangle) for the ground state rotational 
band also agrees very well with theoretical predictions for the negative par- 
ity band. The obtained good agreement supports the interpretation of the 
level scheme of ^^^Ce nucleus presented in Ref. [5]. 
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